Structure based screening and molecular docking with dynamic simulation of natural secondary metabolites to target RNA-dependent RNA polymerase of five different retroviruses

Viral diseases pose a serious global health threat due to their rapid transmission and widespread impact. The RNA-dependent RNA polymerase (RdRp) participates in the synthesis, transcription, and replication of viral RNA in host. The current study investigates the antiviral potential of secondary metabolites particularly those derived from bacteria, fungi, and plants to develop novel medicines. Using a virtual screening approach that combines molecular docking and molecular dynamics (MD) simulations, we aimed to discover compounds with strong interactions with RdRp of five different retroviruses. The top five compounds were selected for each viral RdRp based on their docking scores, binding patterns, molecular interactions, and drug-likeness properties. The molecular docking study uncovered several metabolites with antiviral activity against RdRp. For instance, cytochalasin Z8 had the lowest docking score of –8.9 (kcal/mol) against RdRp of SARS-CoV-2, aspulvinone D (–9.2 kcal/mol) against HIV-1, talaromyolide D (–9.9 kcal/mol) for hepatitis C, aspulvinone D (–9.9 kcal/mol) against Ebola and talaromyolide D also maintained the lowest docking score of –9.2 kcal/mol against RdRp enzyme of dengue virus. These compounds showed remarkable antiviral potential comparable to standard drug (remdesivir –7.4 kcal/mol) approved to target RdRp and possess no significant toxicity. The molecular dynamics simulation confirmed that the best selected ligands were firmly bound to their respective target proteins for a simulation time of 200 ns. The identified lead compounds possess distinctive pharmacological characteristics, making them potential candidates for repurposing as antiviral drugs against SARS-CoV-2. Further experimental evaluation and investigation are recommended to ascertain their efficacy and potential.


Introduction
Everyone has gained a lot of knowledge about viruses after the Covid-19 outbreak over the past three years.The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), human immunodeficiency virus (HIV), hepatitis C virus (HCV), dengue, and Ebola viruses, belong to Retroviridae family, having RNA genome between 7 and 11 kb in size, are notable causative agents of major human diseases [1,2].The enzyme known as reverse transcriptase, which was independently identified in 1971 by American virologists Howard Temin and David Baltimore, gave rise to the term of "retrovirus" [3,4].The dual genetic feature enables the virus to spread from cell to cell as packaged RNA while leaving a copy of the DNA in each infected cell's chromosomes, where it can pass from one generation to the next [5].Research on retroviruses has significantly influenced the fields of genetics, molecular biology, biotechnology, medicine, cellular growth control and carcinogenesis.
The novel SARS-CoV-2 along HIV, HCV, dengue, and Ebola viruses, are currently without a reliable cure.In order to prevent the disease from spreading further, efforts were made too quickly to diagnose the infected individuals and subsequently isolate them [4].Authorized medications and therapies are under clinical trials, and secondary compounds extracted from microbes and medicinal plants are all being explored continuously to cure infections caused by these retroviruses [6].Scientists from various countries are currently striving to develop a vaccine against these viruses, but developing vaccines for RNA viruses, is a very difficult task.
Traditional treatments are mostly composed of natural secondary metabolites which are extracted from plants and microbes [7].They are an important repository of essential chemicals required for the invention of novel pharmaceutical drugs, however, their application has declined in the previous 25 years due to practical challenges in investigating them in molecular target specific experiments [8].These technical hurdles to natural metabolites assessment are now being diminished as new advanced approaches are developed.The therapeutic value of a phytochemical is directly proportionate to the mechanism by which it binds to its cellular target [9].It is the most logical and important step in the medication development process.In retrovirus therapy, there are several pharmacological targets to inhibit, including the viral envelope protein, RNA, reverse transcriptase, integrase, and main protease enzyme [10,11].
A cost-effective, time-efficient, and advantageous practical screening strategy is in silico target identification for metabolites based lead molecules [12].The two most often employed in silico methodologies are quantitative structure activity relationship (QSAR) and molecular docking [13], whereas the molecular docking is more beneficial and fast approach for the screening of target macromolecules (protein, DNA, and RNA) [14].The method also helps to study the mechanisms of action of natural compounds with their target proteins, as well as facilitating the invention of novel derived medications.Ligand-protein docking has resulted in the development of several secondary metabolite-based medications [10,15].Computational tools and software play vital role in drug discovery by accelerating the virtual screening, provide detailed simulation of drug-target interaction and precise molecular modeling [16].These tools and software assist to identify the biological target and their validation, predict pharmacokinetics and pharmacodynamics and helps to optimize the lead compound(s) [17].
The RNA-dependent RNA polymerase (RdRp) participates in the synthesis, transcription, and replication of viral RNA, control the synthesis of RNA using an RNA template.RdRp has been recognized as a possible target for the development of anti-SARS-CoV-2 and other viral therapeutics [18].The antiviral drugs based on nucleotide analogues that target the outer components of the RNA virus's reproduction and functioning machinery [10,19,20].Recently, natural products from marine source were evaluated as an anti-ZIKV agent by targeting RdRp [21].The current study was therefore aimed to inhibit viral transcription, replication, and maturation by targeting RdRp with plant, fungal and bacterial derived secondary metabolites.Furthermore, comparing multiple interactions between ligand-receptor complexes would aid researchers in locating more exact target sites in receptor proteins.The combined targeting of numerous viral proteins with active residues can reveal a single leading treatment candidate against multiple proteins, which can be of immense help in controlling the infections in a timely manner.It will open the way for the investigation of further possible medication candidates with high bioavailability and efficacy against infections caused by retroviruses.

Ligand structure acquisition
Literature and secondary metabolite databases were explored for the collection of compounds with antiviral potential.The 3D structures of 200 antiviral secondary metabolites were retrieved from NCBI's Entrez PubChem database (https://pubchem.ncbi.nlm.nih.gov/;accessed on 07 March 2023) in SDF format.The phytochemicals were energy minimized and then converted into pdbqt format through Open Babel tool [22] of PyRx.The 3D chemical structures of remdesivir (PubChem ID 121304016), abacavir (PubChem ID 441300), and sofosbuvir (PubChem ID 45375808) were used as standards in molecular docking studies.

Target receptor preparation
RNA dependent RNA polymerase (RdRp) of five human pathogenic viral species were selected as target or receptor proteins in this study.Protein Data Bank (PDB) was used to find the molecular structures of RdRp of SARS-CoV-2 (PDB ID: 7B3B) [20], HIV-1 (6UK0) [23], Hepatitis C (4OOW) [24], Ebola (7YER) [25], and dengue virus (5K5M) [26] as receptor proteins.The receptor protein preparation includes adding hydrogen atoms, removing any extra ligand, solvent residues and water molecules, 3D protonation, energy minimized and then converted into pdbqt format through Macromolecule tool of PyRx.Protein preparation was carried out with AutoDock tool Vina [27] of PyRx allocated the rotatable bonds in the ligands.AutoGrid was utilized to create the grid parameters and dimensions.PROCHECK Ramachandran plots were used to assess the quality of protein structures.

Molecular docking
Molecular docking was performed to screen strong antiviral secondary metabolites against RdRp of five different viruses using PyRx virtual screening software [27].After prepration, the protein was imported into PyRx, converted to pdbqt format, and assigned as the receptor protein.A grid box was defined with AutoGrid dimensions of (X: 100. to increase the inquiry into the binding conformations of receptor proteins and ligands.The value of exhaustiveness was set to its default (i.e., 8).The docked poses with the lowest root mean square deviation (RMSD) and binding free energy were chosen.A significant evaluation criterion for screening prospective drugs against RdRp enzyme is the docking score between the receptor protein and ligand [28].Lastly, the visualization of interactions between the receptor protein and important active compounds was done using the Biovia Discovery Studio 2021 [29] software.

Drug scanning
The selected compounds were analyzed by an online tool available at (http://www.scfbio-iitd.res.in/software/drugdesign/lipinski.jsp/;accessed 15 April 2023) to check drug-likeness of these compounds with strong docking scores.Lipinski's rule of five (Ro5) [30] identifies druglike or non-drug-like compounds and predicts the success or failure of drug candidates in the preclinical stage based on five rules.According to Ro5, the molecular mass of the compound should be <500, hydrogen bond donors should be �5, hydrogen bond acceptors should be �10, log P value should be �5, and molecular refractive index should be between 40-130 [30].

Molecular dynamics simulation
One secondary metabolite against each receptor protein with best docking score and druggability properties was selected for molecular dynamics (MD) simulations.The MD simulation analysis of remdesivir (as a control) and RdRp of SARS-CoV-2 complex was also performed and the results were compared with those of best selected compounds.The MD simulation was run for 200 ns using Desmond (Schro ¨dinger LLC) software for MD modeling in replicates [32].First and last pose 2D images of MD trajectories were also obtained for a better comparison.The binding state of the ligand is predicted in static situations with the help of molecular docking and therefore only the earliest phase of the complex is used.The molecular docking is also important as it gives a static view of the bound ligand in the protein active site [33].Newton's classical equation of motion is integrated in MD simulation study which computes the movements of different atoms over time.In the current study, the MD simulation was run to predict the binding status of the ligand-protein complex in the physiological environment [34,35].
The Maestro or Protein Preparation Wizard was used for the optimization and minimization of the ligand-protein complex.The System Builder tool was used to prepare all systems and a solvent model (i.e., transferable intermolecular interaction potential 3 points (TIP3P)) was selected with an orthorhombic box [36,37].In the MD simulation study, the OPLS 2005 force field was employed and the counter ions were introduced for neutralizing the models [38].NaCl (0.15 M) was added for representing the physiological conditions [36].The NPT ensemble with 310 K of temperature and 1 atm of pressure was adjusted for the time of simulation.The models were relaxed before starting simulation and the trajectories were saved after every 100 ps for examination.The protein RMSD values were compared with those of ligands over simulation time to verify the stability of MD simulation.

Screening through molecular docking
For the initial virtual screening, 200 secondary metabolites were docked to RdRp of five selected target proteins.By integrating with active site residues, the protein-ligand docking binding pattern and docking scores filtered out the top five hits as strong active drug candidates against RdRp enzyme.All of the best-chosen secondary metabolites in this investigation demonstrated potent interactions with the active site residues of the catalytic sites of the RdRp enzyme (Table 1).The secondary metabolites cytochalasin Z8 had lowest docking score of -8.9 (kcal/mol) against RdRp of SARS-CoV-2 (7B3B) comparable to standard (remdesivir: -7.The compound talaromyolide D established conventional hydrogen bonds with GlyA:557, AspA:559, and GlyA:449 residues of RdRp enzyme of hepatitis C with a docking score of -9.9 kcal/mol (Table 1, Fig S1).The binding affinity of screened natural metabolites was much lower than all three selected standard chemical drugs approved for the viral treatment by targeting RNA dependent RNA polymerase enzyme (S2

Drug-likeness and ADMET properties
All the selected secondary metabolites followed the drug-likeness properties (Ro5) with no violations except corosolic acid, imatinib, isocochliodinol, and vorapaxar (Table 2).The ADMET based attributes of drug candidates for the biological availability, metabolism, and toxicity were predicted using admetSAR tool in order to further evaluate the pharmacological potential of these secondary metabolites.For a possible medication, a candidate compound that passes the ADMET assessment can be regarded as safe.The findings showed that the selected compounds were non-carcinogens.Overall, the chosen ligands were determined to be suitable as potential medication candidates (Table 3).

Molecular dynamics simulation analysis
Molecular dynamics (MD) simulation analysis was performed to verify the stability of the best ligand against each receptor protein on the basis of findings of molecular docking and druggability analysis.All ligands were found to be firmly bound within the binding pockets of their respective receptor proteins (Fig 7).The ligand-protein complexes did not show drastic fluctuations in their RMSD values throughout simulation period of 200 ns.In ligand-protein complex of RdRp of SARS-CoV-2 and cytochalasin Z8, the protein was found to be stable throughout simulation time while the ligand fluctuated a little after 90 ns from 4Å to 8Å which is quite acceptable.The RdRp of HIV-1 and aspulvinone D complex was revealed to be sufficiently stable throughout the simulation time of 200 ns.The ligand-protein complex of RdRp of hepatitis C and talaromyolide D showed fluctuations in the RMSD values of both protein and ligand till 30 ns but later both gained stability for the rest of the simulation time.Both protein and the ligand in RdRp of Ebola with aspulvinone D complex were overlapped throughout simulation time and the RMSD values showed that the complex was stable throughout simulation time.
In ligand-protein complex of talaromyolide D and RdRp of dengue protein, the RMSD value for Cα-atoms of protein was in the range of 2 to 4Å with a change of 2Å which is acceptable, however, the ligand talaromyolide D showed fluctuations in RMSD values till 82 ns and then the ligand attained stability for the rest of the simulation time.As a control, an antiviral drug remdesivir was also used in the MD simulation study complexed with RdRp of SARS--CoV2.The complex attained stability after 10 ns and the ligand was found to be firmly attached to the receptor protein for the rest of the simulation time.All selected ligand-protein complexes showed stable binding through simulation time and the ligands were firmly bound within the original pockets of their respective target proteins.
The analysis of root mean square fluctuation (RMSF) of receptor proteins complexed with their best selected ligands showed stable peaks throughout simulation time of 200 ns (Fig 8).The findings showed that the ligand-protein complexes are stable.
The PLIF (protein ligand interaction fingerprints) of RdRp of SARS-CoV-2 and cytochalasin Z8 complex before and after simulation revealed that amino acids Asn496, Asn497, Lys500, Ser501, Tyr516, Arg569, and Ala685 are involved in making hydrogen bonds with the ligand and are therefore important to stabilize the complex (Fig 9a).In RdRp of HIV-1 and aspulvinone D complex, the amino acids Arg358, His361, Thr362, Glu370, Tyr405, Gln407, Gly504, Asn418, and Thr419 were found to be involved in stabilizing the complex through Hbond formation ( The complexes of the best selected ligand molecules with five RdRp proteins were also found to be stable through simulation time of 200 ns through the analysis of hydrogen bond interaction stability (Fig 10).The analysis validates the selected ligands to be used as potential candidates against five selected viral proteins.
For the course of trajectory, the stacked bar charts were standardized.According to these stacked bar charts, the contact value of 1.0 shows exact interaction retained for 100% of  The MM-GBSA or molecular mechanics generalized born surface area was also determined for MD simulation (Fig 12).For cytochalasin Z8 and RdRp of SARS-CoV-2 complex, the average ΔG was -49.39 kcal/mol, the ΔG range was -52.91 to -46.89, and the value of standard deviation was 5.42.The average ΔG for aspulvinone D complex with RdRp of HIV-1 protein was found to be -62.50 with its range from -48.19 to -60.94 and standard deviation of 10.67.The complex of talaromyolide D with RdRp of hepatitis C showed average ΔG value of -147.10 with range from -161.34 to -152.86 and standard deviation of 25.18.The average ΔG value for aspulvinone D complexed with RdRp of Ebola virus was found to be -52.73 with range form -50.32 to -41.26 and standard deviation of 7.75.Similarly, for talaromyolide D and RdRp of dengue virus complex, the average ΔG value was -257.41 and the standard deviation value was 17.44 and the range of ΔG was found to be from -298.06 to -229.16.As a control, remdesivir and RdRp of SARS-CoV-2 complex, the average ΔG was -47.77 kcal/mol, the ΔG range was -57.50 to -44.61, and the value of standard deviation was 7.65.The majority of lower binding energy values of best selected ligands with their receptor proteins are validating that the residues of all receptor proteins contributed significantly in stabilizing the ligand molecules near their binding sites.

Discussion
Prior to performing laboratory trials, researchers can analyze the binding mechanisms of different ligands with the aid of novel computational techniques.Modern drug discovery methods such as docking of molecules and dynamic simulation provide essential information on the stability and active site targeting of possible medicinal products against certain receptor proteins [17].These bioinformatics tools help to identify the target molecules and their validation as well as enhance the accuracy and efficiency of treatment to cure these fatal viral diseases in very short time.The secondary metabolites such as alkaloids, phenolics, terpenoids and other bioactive compounds have proven their antiviral activity due to their vast structural efficacy [39].The molecular docking method has been used to demonstrate the ability of several candidates for ligand molecules to bind to various receptor proteins.Various docking strategies have been used depending on the type and degree of analysis as well as the receptor-ligand type.The interactions, energy verifications, and patterns of ligand binding to the receptor moieties' catalytic site are carried out by the docking studies.The testing of various collections of ligand compounds towards various targets, such as proteins, can be aided by computer-assisted drug discovery [40].Maximum binding pocket occupancy, powerful H-bonds, and a minimal energy structure support the possibility of ligand molecules interacting with receptor protein active residues.Before beginning experimental research, computational studies assist scientists in estimating the outcomes of the prospective study.
Diverse geographic areas and the human population have employed wild flora for therapeutic purposes for many years.Natural secondary metabolites used to maintain health and treat a variety of disorders have been found to include a wide spectrum of chemical substances with bioactive components [2].Secondary metabolites with a variety of biological roles are produced by plants and microorganisms.These secondary metabolites fall into categories based on numerous characteristics, including chemical structure (the existence of rings or sugars), composition (the presence or absence of nitrogen), solubility in organic solvents or water, and metabolic process.Carotenoids, phytosterols, phenolics, alkaloids, nitrogen compounds, and organosulfur compounds are among them.Each class is subsequently subdivided into subclasses.Phenolics, in particular, are a diverse category of chemicals, with the phenolic hydroxyl groups serving as the common structural component.These are commonly found conjugated with sugars and organic acids and are classified into seven major groups (i.e., phenolic acids, hydrolysable tannins, flavonoids, stilbenes, coumarins, curcuminoids, and lignans) [9,41].
According to recent research, some of these natural substances can lower the risk of numerous chronic illnesses caused by retroviruses, because of their antibacterial, antioxidant, and anti-inflammatory and immunomodulatory properties [42].Recent breakthroughs in analytical methodologies, such as synthetic biology, metabolic engineering, and metabolomics, as well as computer tools and databases, are providing significant tools for natural chemical drug discovery.Natural products have been shown to have significant antiviral activity, typically against many viral families, implying that they could be beneficial as broad-spectrum antiviral medicines [43].A drug's antiviral activity can influence viral DNA and RNA synthesis, viral reproduction or virus entry and is influenced by viral structure and the replication cycle.The antiviral activity of secondary metabolites can be tested using several biological assays to determine cytotoxicity, cytopathic effect suppression, and the ability to stop viral spread from cell to cell, hence restricting and/or combating viral circulation.Some natural elements present in plant-based foods and dietary supplements show effectiveness against SARS-CoV-2 via affecting cytokine production, modulating cell signaling pathways relevant to inflammation, and even directly interacting with virus targets [44].This has been established through the use of in vitro and in silico molecular docking studies to determine binding strength.Furthermore, the aflavin from Camellia sinensis (tea) has pharmacological activity against SARS-CoV-2 via binding to RNA-dependent RNA polymerase [45,46].In this study, aspulvinone D showed comparable significant interaction with RdRp of SARS-CoV-2.Apigenin has been shown to inhibit SARS-CoV 3CLpro, whereas luteolin adheres to the exterior of the SARS-CoV spike protein and prevents virus entrance into host cells [10].In our previous study [10], plant derived secondary metabolite pedalitin has binding energy of -13.832 kcal/mol against RdRp enzyme of SARS-CoV-2 and formed conventional hydrogen bonds with Tyr619 and Thr556 residues.It has been reported that the quinoxaline derivatives have excellent binding configuration and affinity as well as stability with HIV reverse transcriptase enzyme [47].These results are comparable with the current findings for the binding interaction of aspulvinone D with RdRp stability of HIV-1 and Ebola virus, as well as stability of the complex in dynamic simulation studies.
Various plants and phytochemicals have anti-HIV potential.Gallic acid and ellagic acid have anti-HIV properties and may operate as post-entry inhibitors by decreasing HIV protease activity [48].In vitro, hydroxytyrosol inhibited the viral integrase enzyme and the fusion of the viral envelope with host cells [49].Apigenin, together with epicatechin, (-)-epicatechin-3-Ogallate (ECG), and epigallocatechin-3-gallate (EGCG), have antiviral properties against HIV, likely through blocking the HIV-1 protease enzyme and reverse transcriptase [45,46,49].In addition, studies have examined the antiviral action of several secondary metabolites against the hepatitis C virus (HCV), that is one of the leading causes of chronic liver disease [45,50].Caffeic acid and the flavone apigenin have been shown to impede HCV replication, whilst quercetin-3-O-rutinoside (rutin), curcumin, epigallocatechin-3-gallate (EGCG), and gallic acid have been shown to hinder HCV entrance [51][52][53].Different alkaloids extracted from Atropa belladonna have been studied for the treatment of EBOLA virus.There is no single drug that is universally accepted for treating dengue fever.However, extracts obtained from various plants are expected to reduce the severity of dengue viruses by increasing platelet count, including Carica papaya (papaya) leaf extract used by the local people of India to accelerate the healing process [54,55].Some clinical research has shown that giving C. papaya leaf extract to patients with dengue fever boosts their platelets count [56].
Molecular dynamics simulation is a computational method used to analyze the dynamic interactions within complex systems comprising atoms and molecules over time [57].Structural parameters such as Root Mean Square Deviation (RMSD), Root Mean Square Fluctuation (RMSF), and the number of intermolecular hydrogen bonds are commonly employed to evaluate the stability, dynamic behavior, and compactness of protein-drug complexes.Specifically, RMSD of the protein backbone serves as a key indicator of stability, with lower values indicating a more stable complex.In this study, we utilized RMSD to assess the stability of RdRp in conjunction with specific phytochemicals, comparing it to the RMSD of the free protein.Mean RMSD values for both free RdRp and RdRp in complex with the mentioned compounds were quantified.Our results indicate that the phytochemicals did not significantly perturb the structural stability of these enzymes.Moreover, molecular dynamics simulation study assists as a vital computational instrument to examine the different chemical, physical and biological processes at molecular level.The finding of this analysis is aligned with a previously reported study [39].
Hydrogen bonds play a crucial role in the formation of secondary and tertiary structural protein motifs [58].The binding affinity of a drug towards a protein target in molecular dynamics simulations can be explained by the formation of H-bonds between a ligand and a protein motif, where a higher number of H-bonds indicates stronger interactions.In the case of RdRp, cytochalasin Z8, aspulvinone D, and talaromyolide D exhibited hydrogen bonds with receptor proteins.The high affinities of the ligands to RdRp and the stability of these complexes over time can be attributed to the formation of hydrogen bonds.Small non-covalent or reversible covalent inhibitors, like phytochemicals, are generally preferred due to their advantages in terms of side effects and toxicity when compared to covalent inhibitors.
The current study delivers valuable insights into the antiviral potential of secondary metabolites derived as inhibitors of RNA-dependent RNA polymerase via computational docking and MD simulation methods.Additional validation of the results by in vitro assays is necessary to verify the effectiveness and specificity of these drug candidates against RdRp.Moreover, potential challenges including off-target effects, toxicity, and bioavailability needs to be explored via experimental validation.In-depth in vitro as well as in vivo analysis should be carried out to further validate the in silico findings, and evaluate the safety and therapeutic potential of selected metabolites in future.

Conclusion
The ongoing COVID-19 outbreak, and other viral diseases encounter a tremendous effect on the world economy.This study aimed to identify potential inhibitors against key viral enzyme RNA-dependent RNA polymerase.We screened 200 extracellular secondary metabolites from plants and microbes known for their potency against various viral protein receptors.Top five secondary metabolites were screened out with consideration against the RNA-dependent RNA polymerase of SARS-CoV-2, HIV, hepatitis C, Ebola and dengue virus according to hydrophobic associations, binding amino acid residues and lowest binding affinities.Cytochalasin Z8 had the lowest docking score of -8.9 (kcal/mol) against RdRp of SARS-CoV-2 interacted with the binding pocket residues.In this study, aspulvinone D and talaromyolide D demonstrated broad-spectrum antiviral potential against multiple viruses.Aspulvinone D showed the lowest binding affinity with the RNA-dependent RNA polymerase (RdRp) of HIV-1 and Ebola virus, with docking scores of -9.2 and -9.9 kcal/mol, respectively.Talaromyolide D showed interaction with RdRp of hepatitis C and dengue virus, with docking scores of -9.9 and -9.2 kcal/mol, respectively.Drug-behavior and ADMET profiling displayed that selected natural secondary metabolites were found to be non-toxic and non-carcinogenic and possess strong antiviral effects against viral RNA-dependent RNA polymerase enzyme.Consequently, these compounds may be taken into account for further experimental trials and investigations to explore their antiviral potential to treat viral infections in future.

Fig 7 .
Fig 7. Root-mean-square deviation (RMSD) values of selected receptor proteins with their best ligands.(a) RMSD of the C-alpha atoms of RdRp of SARS-CoV-2 and cytochalasin Z8; (b) RMSD of the C-alpha atoms of RdRp of HIV-1 and aspulvinone D with time; (c) RMSD of the C-alpha atoms of RdRp of hepatitis C with talaromyolide D; (d) RMSD of the C-alpha atoms of RdRp of Ebola with aspulvinone D; (e) RMSD of the C-alpha atoms of RdRp of dengue with talaromyolide D; (f) RMSD of C-alpha atoms of RdRp of SARS-CoV2 and remdesivir as a control.The variation of protein RMSD is shown on the left Y-axis through time.The variation of ligand RMSD is shown on the right Y-axis through time.https://doi.org/10.1371/journal.pone.0307615.g007